stdev <- function(data, na.rm = TRUE) {
	return (sqrt(variance(data, na.rm)))
}